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Abstract 

We describe the program package gentle/4f an, which consists of the Fortran codes gentle_4f an.f , 
4fan.f, and gentle_tic_qed.f and is devoted to the description of e + e~ annihilation into four 
fermions. The codes are based on the semi-analytical approach. Initial state QED corrections are 
taken into account. The program versions, which are described here were used in the 1995 workshop 
'Physics at LEP 2'. 

> Contributions to the working groups on 'Event Generators for WW Physics' and 'Event Gen- 
erators for Discovery Physics' of the 1995 Workshop on Physics at LEP 2, to appear in the 
proceedings (G. Altarelli, T. Sjostrand and F. Zwirner (eds.), CERN 96-01) 



1 Supported by European Union with grant INTAS-93-744. 



1 Introduction 



This note collects two slightly updated contributions to the working groups on 'Event Genera- 
tors for WW Physics' jl] and on 'Event Generators for Discovery Physics' of the Workshop 
on 'Physics at LEP 2' held at CERN in 1995. The underlying formulae may be found in refer- 
ences |§-[P]. In this series of papers we advocate a semi-analytical approach to the reaction 

e + e" /i/i/ 2 / 2 (n 7 ). (1) 

The approach is characterized by the attempt to perform at least part of the eight-dimensional 
phase space integration analytically. As a result, the numerical integrations are smoothened 
and often faster than in a completely numerical approach. Further, the numerical accuracy 
is well under control. The analytical integrations are performed over (part of) the angular 
variables, while the boson invariant masses are left for numerical integration. 

Four-fermion production processes may be very complex and proceed via a large variety of 
intermediate states. Doubly resonant processes are usually calles 'basic processes'. The others 
are called background contributions. We repeat the classification of reference ||, with the 
number of Feynman diagrams in the charged current (CC) classes are shown in table |I| and for 
the final states corresponding to the neutral current (NC) classes in table 0. 
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Table 1: Number of Feynman diagrams for class CC final states. 
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Table 2: Number of Feynman diagrams for class NC final states. 
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Description of the package 

The GENTLE/4f an package is designed to compute selected total four-fermion production cross- 
sections and final-state fermion pair invariant mass distributions for charged current (CC) 
and neutral current {NC) mediated processes within the Standard Model (SM). For the CC03 
subprocess, the W production angular distribution is also accessible. In the NC case, SM 
Higgs Production is included. The phase space integration is carried out by a semi-analytical 
technique, which is described below. The GENTLE/4fan package is written in Fortran. It 
consists of two branches. The basic branch gentle_4f an.f contains all features of the package 
but complete initial-state radiation (ISR) to NC processes. The subroutine FOURFAN, which 
is called by gentle_4f an.f performs the computation of NC cross-sections and is described 
in the next section. The independent branch gentle_nc_qed.f includes complete ISR to NC02 
and NC08 and will be merged into gentle_4f an.f . 

Program features: 

1. Method of integration: 

The package is a semi-analytical one. Without (with) ISR, the phase space is parame- 
trized by five (seven) angular variables and the final-state fermion pair invariant masses 
(plus the reduced center of mass energy squared). All angular phase space variables 
are integrated analytically. The resulting formulae are input to the package. Invariant 
masses are subsequently integrated numerically with a self-adaptive Simpson algorithm. 
Optionally, for the CC03 subprocess, the W production angle may also be numerically 
integrated. The method is numerically stable and, in most cases, very fast. 

2. Possible final states: 

The package may treat all four-fermion final states which do not contain identical particles, 
electrons, or electron neutrinos. This means that the package accesses all final states that 
are described by annihilation and conversion type Feynman diagrams: 
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(1) CC03 (with complete ISR) § 

(2) NC02, NC08 (with complete ISR) 

(3) CC9, CC10, CC11 @ 

(4) NC06, NC10, NC24, NC32 § 

(5) NC + Higgs [0 

Via flags, cross-sections for subsets of Feynman diagrams may be extracted. 

3. Cuts 

Cuts may be imposed on invariant masses of fermion pairs and on the invariant mass of the 
final-state four-fermion system. Using the structure function approach in gentle_4f an . f , 
cuts on the electron/positron momentum fraction can be imposed. For the CC03 sub- 
process, cuts on the W production angle are enabled. 

4. Initial state radiation 

ISR is implemented into the package. Universal ISR is present for all processes [jSJ. In 
addition, the package includes complete, i.e. universal and non-universal ISR for the 
CC03, NC02, and NC08 processes [[| [|. Non-universal ISR does not contribute to 
annihilation diagrams. It may be argued that non-universal ISR is very small, O(10~ 3 ), 
for conversion- annihilation interferences. The speed of the package is reduced, if non- 
universal ISR is included, due to its complex analytical structure. 

5. Final state radiation 

Final state radiation is not implemented. 

6. Treatment of final state decays 

Final state decays are not accounted for. 

7. Treatment of the Coulomb Singularity 

The Coulomb singularity is included according to reference ||]. 

8. Treatment of the Anomalous Couplings 

Anomalous couplings are not included. 

9. Treatment of masses 

In general, final-state masses are neglected in the matrix elements. Where needed, how- 
ever, masses are retained in the phase space. In addition, masses of heavy particles 
coupling to the Higgs boson are taken into account where appropriate; see section 3. 

10. Hadronization 

No interface to hadronization is foreseen. 
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Input parameters 



All input parameters are set inside the Fortran code. gentle_4f an.f uses the following flags, 
set in the subroutine WWINOO: 

IBCKGR: CC03 case (IBCKGR=0) or CC11 case (IBCKGR=1) 

IBORNF: Tree level (IB0RNF=0) or ISR corrected (IB0RNF=1) quantities 

ICHNNL: CC03 (ICHNNL=0), CC11 with specific final state [^^(ICHNNL = l),lvqq 

(ICHNNL = 2,3), (ICHNNL = 4)], and inclusive CC11 (ICHNNL=5) 

ICOLMB: Inclusion of Coulomb singularity (IC0LMB=1,...,5) or not (IC0LMB=0) 

Recommended value: IC0LMB=2 
ICONVL: Flux function (IC0NVL=0) or structure function approach (IC0NVL=1) for 

ISR. Recommended value: IC0NVL=0 
IGAMZS: Constant Z width (IGAMZS=0) or s-dependent Z width (IGAMZS=1) 
IINPT: Input for tuned comparison [|l(| (IINPT=0) or preferred input (IINPT=1) 
IIQCD: Naive inclusive QCD corrections are included (IIQCD=1) or not (IIQCD=0) 
IMMIM: Minimum number of a moment requested by IREGIM 
IMMAX: Maximum number of a moment requested by IREGIM 
IONSHL: On-shell (I0NSHL=0) or off-shell heavy bosons (I0NSHL=1) 
IPROC : CC case (IPR0C=1) or NC case (IPR0C=2, call to FOURFAN is initialized) 
IQEDHS: Determination of the universal ISR radiator: 
0(a) exponentiated (IQEDHS=-1,0); 

0(a) exponentiated plus different 0(a 2 ) contributions (IQEDHS=1,...,4) 
Recommended value: IQEDHS=3 

IREGIM: Calculation of the total cross-section (IREGIM=0), the moments of the radi- 
ative loss of final-state four-fermion invariant mass (IREGIM=1), the moments 
of the radiative energy loss (IREGIM=2), the moments of the W mass shift 
{V^+~^ ~ 2-Ww) (IREGIM=3), and the first moments of cos (n6 w ), 
n = 1, ...,4 (IREGIM=4) 

IRMAX : Maximum value of IREGIM 

IRSTP : Step in a DO loop over IREGIM 

ITVIRT: Non-universal virtual ISR included (ITVIRT=1) or not (ITVIRT=0) 
ITBREM: Non- universal bremsstrahlung included (ITBREM=1) or not (ITBREM=0) 
IZERO : See equation (4.5) of ||. Recommended value: IZER0=1 
IZETTA: See equation (4.21) of ||. Recommended value: IZETTA=1 

In the gentle_nc_qed.f branch, only the flags IBORNF, IONSHL, ITVIRT, ITBREM are used. 
The additional flag IBOSON in gentle_nc_qed. f distinguishes between the NC02 and the NC08 
processes. 

The center of mass energy squared is chosen by setting the variable IREG and the parameters 
ISMAXA or ISMAXB in the main program. The following input may be changed by the user: 
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GFER 
ALPW 
AME 
AMZ 
AMW 
GAMZ 
ALPHS 



G, 
a(2M w ) 
m e 
M z 
M w 

T z 

a s (2M w ) 



1.16639 xlO 5 GeV 2 , the Fermi coupling constant 

1/128.07, the running fine structure constant at 2M W 

0.51099906 xl0~ 3 GeV, the electron mass 

91.1888 GeV, the Z mass, 

80.230 GeV, the W mass 

2.4974 GeV, the Z width 

0.12 



Output 

The following derived quantities are computed in gentle_4f an.f and printed in the output: 
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GVE and GAE are the electron vector and axial vector couplings, GWF is the fermion-M^ coupling, 
and IGWWGI and |GWWZ| are the trilinear gauge boson couplings for the photon and the Z 
respectively. Further the output repeats the flag settings. After the cross-section calculation, 
the following output is printed: 



SQS = 

XSEC0 = cr tot (s) in nanobarns 



(2) 



In addition, the calculated MOMENTS are printed. In the first column IREGIM is printed. The 
second column is arranged in blocks of three lines each. The first line contains the integer n. 
The second line contains the n th moment of the physical quantity indicated by IREGIM. The 
third line contains the dimensionless n th moment obtained through division of the n th moment 
by the proper power of \fsj2. 

Although variable names are slightly different, gentle_nc_qed. f uses the same derived quan- 
tities as gentle_4f an.f . For one run, gentle_nc_qed.f outputs the used flag values together 
with the fermion code numbers I FERM 1 / 1 FERM2 , the color factors RNC0U1/RNC0U2, the masses 
AM1/AM2, and the invariant pair mass cuts CUTM12,CUTM34 for the final-state fermion pairs. In 
addition, the lower cut CUTXPR on the ratio of the four-fermion invariant mass squared over the 
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center of mass energy squared, s'/s is output. The main output, however, is an array of center 
of mass energies and the corresponding total cross-sections. 

Availability 

The codes and this description are available from the authors upon E-Mail request or via WWW 



gentle_4f an.f 
gentle_nc_qed . f 

DESY-Zeuthen 96-05 



from |http : // www . if h . de/| bardin/gentle_4f an.uu 
from |http : //www, if h. de/| ~ lehner/gentle_nc_qed.uu 



from |http : //www, if h. de/theory/publist . html 
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3 4fan version 1.3 



Authors: 

D. Bardin bardin@ifh.de 

A. Leike leike@graviton.hep.physik.uni-muenchen.de 

T. Riemann riemann@ifh.de 

Description of the package 

4f an is a semi-analytical program which calculates the process 

e + e" - /1/1/2/2, (3) 

where the three involved fermions e, f\ and f 2 must be in different electroweak multiplets (the 
NC32 class) ||; see table 2. Optionally, Standard Model Higgs production can be included 0. 
For calculations at the Born level, 4f an can be used as a stand-alone program. For the calcula- 
tion of cross-sections including initial state radiation, the initial state radiation environment of 
the code gentle_4f an. f must be used which calls FOURFAN as a subroutine. For the description 
of gentle_4f an.f we refer to section 2. In the following, we describe the stand-alone program 
4fan.f: 

Six of the eight integrations of the four particle phase space were done analytically. The 
two remaining integrations over s\ = [p(fi) + vi.fi)} 2 an d s 2 = [p(/2) + vih)] 2 are performed 
numerically allowing the inclusion of cuts for these variables. 

Finite mass effects are taken into account using the following approximations: 



• The phase space is treated exactly. 

• In the Higgs contributions and the conversion diagrams e + e~ — > (77) — > /1/1/2/2, the 
masses are treated up to order 0[m 2 (fi) / Sj\. 

• Fermion masses are treated identically in traces and Higgs couplings. 

• The Higgs width is calculated including the decays into b-, c- and r-pairs. 



The numbers quoted in the tables of the report [|1(J are produced for zero fermion masses except 
in the Higgs couplings. The Higgs propagator is always connected with s 2 by convention. 

The initialization routine BBMMIN contains the input from the Particle Data Group [IT]. In 
the subroutine DSDSHSZ, the interferences between the three main subsets of the diagrams of 
the NC32 class are calculated as well as those with the Higgs signal diagram. Their sum yields 
the double differential cross-section. Single interferences between these subsets are not printed. 

The numerical integration is done by a twofold application of a one-dimensional self-adaptive 
Simpson algorithm with control over the relative and the absolute error. The singularities due 
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to resonant vector boson propagators are eliminated by transformations of the integration 
variables. To avoid numerical instabilities, the kinematical functions resulting from the six-fold 
analytical integration are replaced by Taylor expansions near the borders of the phase space. 
The shortest calculation time is achieved by a choice of the required absolute and relative errors 
in such a way that they give approximately equal contributions to the error of the output. 

The calculation time of a Born cross-section is several seconds on an HP workstation, de- 
pending on the required accuracy and the cuts on s\ and S2- Ten times higher accuracy needs 
approximately two times longer calculational time. 

Input and output are transferred through the arguments of the subroutine only. 

Usage of the program: 
CALL FOURFAN (EPS , ABSE , IF1 , IF2 , S , S1MIN , S 1MAX , S2MIN , S2MAX , AMH , IOUT , OUT) 

Input: 



EPS, ABSE: 
IF1.IF2: 



S: 

S1MIN,S1MAX: 

S2MIN,S2MAX: 

AMH: 

IOUT: 



I0UT=1 : 
I0UT=2: 
I0UT=11, 12: 



The required relative and absolute error. If one of 

the two criteria is fulfilled, the calculation stops. 

Integers specifying the two final fermion pairs as in 

the Monte Carlo particle numbering scheme, see 

Particle Data Group |TTJ, chapter 32. 

The cm. energy squared of the e + e~ pair. 

The integration bounds of Si. 

The integration bounds of S2- 

The Higgs mass. 

Integer, selecting the output. 

Currently I0UT=1, 2, 11 and 12 are implemented: 
Total cross-section a t without Higgs. 
Differential cross-section da/ds2 without Higgs. 
The same as I0UT=1, 2 but with Higgs. 



The units of the input (if required) are GeV 2 or GeV. 
Output: OUT Depends on the value of IOUT. The output is given in fb or in fb/GeV. 



On HP workstations 4f an must be compiled with the -K option. 
Availability 

The code and this description are available from the authors upon E-Mail request or via FTP 
or WWW 



from [ftp : / /gluon . hep . physik . uni-muenchen . de 



4fanvl3.f 

4fanvl3.f from |http : //www, if h. de/theory/publist . html 

DESY-Zeuthen 96-05 from |http : //www, if h. de/theory/publist . html 
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4 Concluding remarks 



Not all of the final-state topologies of tables 1 and 2 have been treated by our approach so 
far. Those which have been treated belong to the classes CC11 and NC32 and are printed in 
boldface in the tables. 

Presently we are studying the semi-leptonic and leptonic processes of the classes CC20, 
NC48, NC21, which are printed in the tables in roman [T2| . 



For a complete treatment of the NC32 class, the four jet production from e + e annihilation 
into qqgg has to be covered in addition to 4/ production ||13||. 



Besides a verification of the Standard Model predictions for WW and ZZ pair production 
and in Higgs boson searches via ZH production, one is also interested in searches for anomalous 
triple gauge boson couplings. A description of CC03, CC11 cross-sections with anomalous 
couplings is under development ]14| . 



For a detailed study of the underlying physics, angular distributions are extremely helpful. 
Although in the semi-analytical approach one has a limited flexibility concerning distributions 
and cuts, we made several attempts to calculate some of them [||, [L4], [I5[| . 
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